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Abstract 

m 

<^ • Motivated by the work of Mersini, the particle production related to the tunneling in false 

vacuum decay is carefully investigated in the thin-wall approximation. It is shown that in this 
case the particle production is exponentially suppressed even when the momentum is comparable 
to the curvature scale of the bubble. The number of created particles is ultraviolet finite. 
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1 General discussion 



^ I False vacuum decay by the tunneling in the WKB approximation is described by the 0(4)-symmetric 
(N) ' bounce solution, <frb{Vx 2 + t 2 )• In the limit of small energy-density difference, e, between the true 
^ ■ and false vacuum the bounce looks like a large four-dimensional spherical bubble of true vacuum 
{Sj . separated by a thin wall from a sea of false vacuum. The evolution of the system after tunneling 
O is described by the 0(3, 1) invariant solution that is obtained by the analytic continuation of the 
<3 bounce to the Minkowskian time r — > it, see 0. Since the field 4> describing the bubble nucleation end 
expansion is coupled to itself and to the fluctuation field, the false vacuum decay through the barrier 
penetration is accompanied by the particle production. Using the formalism developed in [2HH| for the 
Ph. calculation of particle number the particle production by the thin- walled bubble was considered in [I]. 
<— | However, we fail to corroborate the main conclusion of |lj that for the momentum comparable to the 
^ \ curvature scale of the bubble there is a strong enhancement of the particle production. In this paper 
we would like to point out that, in general, the particle production in the thin- wall approximation is 
strongly suppressed. The calculations have been undertaken to exhibit the exponential suppression 
of the particle production, even when the momentum is comparable to the curvature scale of the 
bubble, rather than for actual use in calculation of particle number at all momenta. Throughout this 
paper the metric signature in Minkowski space-time is (— , +, +, +) and units h — c — 1 are used. 
The calculations are present in sec. 2. 



2 Particle spectrum 

Throughout this paper we restrict ourselves to the consideration of thin-wall case. This approxi- 
mation for the bubble formation is considered in pQ. We consider a scalar field (f> defined by the 
Lagrange density 

c = -\{d»<p) 2 -u{<p). (i) 

The potential has the form 

U{(j>) = ^0 2 (0 - 2a) 2 - e(0 3 /2a 3 - 30 4 /16a 4 ), (2) 

where A, a and e are positive parameters. The potential (0) has a local minimum at </>/ = 0, U(0) = 0, 
and a global minimum at <p t = 2a, U(2a) = —e. Bounce solution with a good accuracy is 0& = 



a{l — tanh(/i[p — R])} where p, = a\/A and R = 2/i 3 /eA. Thickness of the bubble wall is L = 2//i. For 
the potential (j2J) the condition for the validity of the thin-wall approximation is R/L = p 4 /e\ ^> 1. 
In other words, in this limiting case the energy- density difference between the true and false vacuum 
is much smaller than the height of the barrier U(<j)). For the calculation of particle spectrum we use 
the standard technique developed in (21 E] . The results relevant for our discussion can be summarized 
as follows. The coordinate system used in the Euclidean region is (p, \i <p) where (6, ip) are usual 
angle coordinates on two-dimensional sphere and (p, x) are related to r = \x\ and r as follows 



r = psin(x), r = -pcos(x), 
< X < 7r/2, < p < oo. 

The coordinate system in the Minkowski region may obtained by replacement (p, x) 
which yields 

r = p M sinh(xM), * = Pai cosIi(xm), 
< Xm < oo, < p M < oo. 

The basic equation which governs the fluctuation field $(£, x) reads 
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(4) 



(5) 



where Dm denotes the partial derivative with respect to pu and L is the Laplacian operator on 
three-dimensional unit hyperboloid. Since the bounce solution has 0(3, 1) symmetry in Minkowski 
region one can label $ by the eigenvalues of the angular momentum, (I, m) as well as by the radial 
momentum, p, and write 

v ( a \^m(Pm) /„x 

Ypt m (x,0,<P) 372—, (6) 

Pm 



where Y v \ m are eigenfunctions of the operator — 1? with eigenvalues 1 +p 2 . After substituting Eq.( 
into Eq.(j3J) and analytically continuing of $ to the Euclidean region by replacement pu — > ip, Xm 
—ix one obtains the equation for the fluctuation field in the under-barrier region 



(P 2 + 1/4) 



U"( 



^ = 0, 



where ip(p) = ^m^p)- This is a zero-energy Schrodinger equation 

K - V(pM = 

of a particle in the potential 

(P 2 + 1/4) 



V(p) 



+ 2/T{3tanlr(/![p - R}) - 1} + e-term. 



(7) 



(8) 



(9) 



Thus, iPm(pm) is given by solving Eq.(jHJ) with p replaced by — ip. Since in the thin- wall approximation 
4>b (0) is very close to the true vacuum value 2a the solution obtained in this way after the analytical 
continuation at the turning point p = will have the form 

Vm(pm) = c l ^)^e^' 2 H^{up M ) + c 2 {p)^e^' 2 Hll ) {up M ), 



where u 2 = U"(4>b(0)). At t — — oo the field is in false vacuum state and correspondingly the 
fluctuation field $ satisfies the vanishing boundary condition when p —>■ oo. The spectrum of created 
particles n{p) is given by [3] 

n[p) = (M P )/e 2 W-ir (10) 

For detailed description of the general formalism see [21 E|- 

If p is of order or smaller than pR one can approximate the potential V as follows 

r v - v/ P 2 P+ < P , 

V = < v 3 (x + v 2 /v 3 ) 2 - k \x\<L, (11) 
[ Ap 2 -v/p 2 p < p_, 

where v = p 2 + 1/4, x = p — R, Vi — v/R 2 + 2p 2 , v 2 = v/R 3 , v 3 = 6/i 4 — 3v/R A , v 3 = v 3 (6 + 
(p/ pR) 2 )/24, k = v 3 (vi/v 3 + v 2 /v 2 ), p± = R ± L. We have neglected the contribution coming from 
the e-term because at = this term equals zero and at <p ~ a this is of order e\/p 2 <C p 2 . The 
general solution of Eq.fjHJ) in regions \p — R\ > L is written in terms of modified Bessel functions. For 
solving this equation we have used the program package Maple 7. 
In region p < p- the solution reads 

ip = a ly /pl ip (2pp) + a 2 ^/pK ip (2pp). (12) 



In region \x\ < L, the exact solution to this equation is expressed in terms of Whittaker W, M 
functions 0. Using the program package Maple 7 one finds that in region \x\ < L the solution has 
the form 

t/j = a 3 W(b, 1/4, z 2 )/J\z] + a 4 M(b, 1/4, z 2 )/ <J\z\, (13) 

where z = v^ 4 (x + v 2 /v 3 ), b = k/Av^ 2 . For s = p/pR by taking into account that in the thin- wall 
approximation (pR) 2 <C (pR) 4 one gets 



k/v^ 



v + 2p 2 R 2 + v 2 /(6p 4 R 4 - 3v))NQp 4 R A - 3v 



and 



((2 + s 2 )p 2 R 2 + v 2 /6p 4 R 4 )/jQp A R A « (2 + s 2 )/V6 



{±2/p + v 2 /v 3 )vl 14 = {±2/pR + v/{6p 4 R 4 - 3v)){6p 4 R 4 - 3^) 1/4 « 

(±2/pR + v/6p 4 R 4 )pRQ 1/4 « ±2 x 6 1/4 . 



It is convenient to introduce the following notations 



W±(s) = W / (&,l/4,^ 2 )/^/|z| 
M±(s) = M(6, \l^z 2 )jJ\z 



±z 



±z 



W'Js) 



(W(b,l/4,z 2 )/ 



±z 



M' ± (s)= \M(b,l/A,z 2 )/J\z 



where ±z = z(x = ±L) » ±0.1304 x (6 + s 2 ) 1 / 4 and the prime denotes differentiation with respect 
to z. These quantities are needed for the matchings at p±. Using the program package Maple 7 one 
can estimate W±, W±, M±, M' ± at the desired value of s. 

In region p > p + the solution satisfying the vanishing boundary condition when p — > oo reads 



i> = a 5y /pK ip (2pp). 



(14) 



For our purposes the quantity of immediate interest is n(p) = |ai(p)/a 2 (p)vr| 2 , see j3]. So, one can 
omit the common multiplier of ai, a 2 coefficients as well as the common multiplier of a 3 , a 4 coefficients. 
The matching at p + gives 



a 3 ~ M' + K ip {2w + ) - M + {K ip (2 fJ ,p + )/2vl / * P+ + K^pp+^p/^ 1 / 4 }, 
a 4 ~ ^ + {^ p (2pp + )/2S 3 1/4 p + + K' ip {2pp + )2p/vl ,A } - W' + K ip {2 W+ ), 



(15) 



the proportionality coefficients are the same for a 3 and a 4 , here and below the prime denotes differ- 
entiation with respect to the whole argument. The matching at p_ gives 

O! ~ {K vp {2pp^/2v l 3 /4 p_ + K' ip (2^)2p/dl /4 }{a 3 W. + a 4 M_} - 

K ip (2pp„){a 3 W'_ + a A M'_}, 

a 2 ~ -{I lv (2pp_)l2vTp- + r ip {2pp_)2p/vl' A }{a 3 W. + a 4 M_} + (16) 

I ip (2pp_){a 3 W'_ + a A M'_}, 

the proportionality coefficients are the same for a\ and a 2 . In the thin- wall approximation pp± ^> 1 
and for evaluating the modified Bessel functions of the 2pp±'s one can use the asymptotic formulae. 
For p <C \J pR one can use the following asymptotic expansions , see jHj, 



I ip (x) ~ exp(x){l + (V + l)/8x + ...}/V2^x, 
I' ip (x) ~ exp(x){l - (3 - Ap 2 )/8x + ...}/y/2^, 

K ip (x) ~ exp(-x){l - (V - l)/8z + ...}^tt/2x, ^ 
^(x) ~ -exp(-a;){l + (3 - 4p 2 )/8x + ...}y/n/2x. 

For p ~ pi? the asymptotic expansions of modified Bessel functions of purely imaginary order and 
their derivatives are derived in [Zj. Motivated in part by the fact that Ki P (x), where x is real variable, 
is the imaginary part of Ii P (x) up to a multiplicative factor, recent investigations have been carried 
out on the real part Ii P (x), denoted by L ip (x), 

L ip (x) = -{h P {x) + I- ip (x)}. 
Any of Li P (x), K ip (x) and Ii p (x) can be constructed from the remaining two functions by the identity 

I ip {x) = L ip {x) - i SmHpn) K ip (x). (18) 

7T 

Both K ip (x) and L ip (x) are real valued and even functions of p. 

For p < x, introducing p = xsin^, where < 9 < tt/2, the leading terms of asymptotic expansions 
are [7j 



Lip(p esc 9) ~ exp(p [cot# + 6'])/ v /2pTcot^ + exp(p [ir - 9 - cot 6»])3p cot 2 9/2n , 
L- p (p esc 6>) ~ exp(p [cot 9 + 9})Jcot 9/2pn + exp(p [it — 9 — cot 9]) tan# sin 9/3n . 



Ki p (p esc 9) ~ exp(— p [cot # + #])y 7r/2pcot , 
i^ p (p esc #) ~ — exp(— p [cot^ + #]) sin^y^Tr cot9/2p. 

For p x the leading terms of asymptotic expansions take the form [7j 

L ip (p) ~ exp(p7r/2){^(0)/2 2 / 3 p 1 / 3 + ^ / (0)/2 1 / 3 p 5 / 3 70} , 
L'iplp) ~ exp(p7r/2){^ / (0)/2 1 / 3 p 2 / 3 - ^(0)/2 2 / 3 p 4 / 3 5} , 
#i P (p) ~ 7rexp(-p7r/2){2 1 / 3 Ai(0)/p 1/3 + 2 2 / 3 Az'(0)/70p 5/3 } 
i^(p) ~ 7rexp(-p7r/2){2 2 / 3 ^ , (0)/p 2/3 - 2 1 / 3 ^(0)/5p 4 / 3 } , 



(19) 



(20) 



where Ai, Bi are Airy functions. 

For p > x, introducing p = x cosh v and a = p(tanh v — v) + 7r/4, where v > 0, the leading terms of 
asymptotic expansions are [7j 



Li P (p sechz/) ~ exp (p 7r/2 ) \/ coth z/{ cos a + sin a }/ y/2p~K , 

L^,(p sechz/) ~ exp(p7r/2) cosh z/{sin a — cosa}/y / 2pn coth z/ , . . 

K ip (p sechz/) ~ exp(— p7r/2)v / 27rcothi/{cosa — sinaj/y/p, 

sechz/) ~ — exp(— p7r/2)v^27rcoshz/{sina; + cos a}/Vp coth z/ . 

We can now analyze the particle spectrum. Using the leading asymptotic terms from Eq.()17|) for 
small values of p one gets 

n{p) « l6(pR) 2 exp(-8pR). (22) 

Notice that this result is different from that one obtained in 0] for small values of p. Taking into 
account that W±, M±, W'±, M' ± are real valued functions from Eqs. (jl5ll6ll8j) one obtains 

_ , » .sinh(p7r) 
a 2 ~ rLe(a 2 ) + % ai, 

7T 

and correspondingly 

(7iKe{a2)/ a.\) + smn (p7r) 
Thus, If p ^> 1 from Eq. (j2^j) one simply concludes that n(p) is suppressed at least as 

n (p) ~ exp(— 2pn). (24) 

As it is clearly seen from Eq. (j23|) n(p) vanishes for that values of p at which Re(a2)/ai becomes 
infinity. Using Eqs. ()19l20l21|) one can easily estimate the term Re(a2)/ai with exponential accuracy. 
If p is large, of order fiR, but small than 2fiR such that arcsin(p/x) is not very close to n/2 from 
Eq. fjTTH) one gets 

Re{a 2 )/a l ~ exp(4p[cot 9 + 9]), (25) 

where < 9 < n/2. For the values of 9 very close to n/2 one has to use Eq. (j20J) . If p is close to 2fiR 
using Eq. (j2~U|) one gets 

Re(d2)/ai ~ exp(47r/i_R). (26) 
For p of order fiR such that p > 2fiR from Eq.()21|) one gets 

Re(a 2 )/ai ~ exp(2p7r). (27) 



For very large momentum, assuming for instance p ^> (fiR) 2 , one can approximate V(p) as follows 

f 4/i 2 - v/p 2 P+ < p, 
V = I -2p 2 - v/p 2 \x\ < L, (28) 
[ 4/x 2 — v/p 2 p < p-. 

Namely, we have replaced the well of the second term of Eq.Q in region \x\ < L by the square 
well, and if p is very large one can neglect the contribution coming from such a deformation of the 
potential in region \x\ < L in comparison with the term v/p 2 . Now in region \x\ < L the general 
solution of Eq. (JHJ) is expressed in terms of Bessel functions of purely imaginary order. With no loss 
of generality one can choose the solution in this region to be real valued function. For instance in 
this region one can take 

ip = a 3 ^pRe{J ip {V2pp)) + a 4y /pRe{Y ip {V2pp)). (29) 



It is clear that for the particle spectrum one gets again the Eq.(|23|) and simply concludes that for 
very large momentum n(p) is suppressed at least as 

n ip) ~ exp(— 2pit). (30) 

The asymptotic behavior of Bessel functions of large purely imaginary order can be found in [H]. 
Namely, in jHj one can find the asymptotic expansions of Ji P (x) and H^\x), 

Jip(x) ~ exp(j)7r/2) exp(? v / P 2 + x 2 — ip arsh(j»/x) — in / 4) / y/2n (p 2 — x 2 ) 1 / 4 , 
H{p\x) ~ v / 2exp(p7r/2) exp(iy / p 2 + x 2 — ipaish(p/x) — m / '4)/7r(p 2 + x 2 ) 1 ^ 4 , 

and using the standard relation Y ip (x) = i{Ji P (x) — H^\x)} can determine the asymptotic behavior 
of Y ip (x). Using Eqs. ()21l31|) one can easily estimate the ratio Re(a 2 )/ai with exponential accuracy. 
The matchings at p + , p_ give 03,4 ~ exp(— pn) and a± ~ exp(— pir/2), a 2 ~ exp(p7r/2) respectively 
Correspondingly, for very large momentum 

Re(a2)/ai ~ exp(2p7r). (32) 

Thus, we have explicitly shown that, in general, the particle production in the thin-wall approxima- 
tion is exponentially suppressed. Using Eqs. (j21l31|) one can also evaluate the ultraviolet finiteness of 
the number of created particles. Without carrying out the explicit computation from these asymp- 
totic formulae one concludes that for very large momentum the particle production is controlled by 
the term n(p) ~ exp(— 2pir) x (rational function of p). As an unessential to this problem we do 

not take into account the terms enclosed in the braces in Eq. (j21)) . For this expression the integral 

00 

/ n(p) is certainly convergent. So we infer that the number of created particles is ultraviolet 

p»(a*R) 2 

finite. This completes the review. 
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